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ABSTRACT 


Aim and Objective: An analysis is made to study the problem of boundary layer forced convective flow and heat transfer of an 
incompressible fluid past a plate embedded in a Darcy-Forchheimer porous medium in presence of transverse magnetic field. 
Thermal radiation term is considered in the energy equation. 

Methodology: The governing partial differential equations are transformed into self-similar ordinary differential equations using 
similarity transformations, which are then solved using Runge-Kutta fourth order method along with shooting technique. 

Results: The numerical results are plotted in some figures and the variations in physical characteristics of the flow dynamics and 
heat transfer for several parameters involved in the equations are discussed. It is found that these parameters have significantly 
effects on the flow and heat transfer. 

Conclusion: It is found that the influence of the physical parameters viz. Porous medium parameter, the Inertial parameter, the 
Magnetic parameter, the Prandtl number, and the Radiation parameter have significantly effects on the flow and heat transfer 
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NOMENCLATURE 

M 

Magnetic parameter 


Constant applied magnetic field 

N 

Radiation parameter 

C P 

Specific heat at constant pressure 

Pr 

Prandtl number 

f 

Dimensionless stream function 


Radioactive heat flux 

k 

Darcy permeability of the porous medium 

s 

Suction parameter. 

k' 

Forchheimer resistance factor 

T 

Temperature of the fluid 


Parameter of the porous medium 

u, V 

Velocity component of the fluid along the x 

k 2 

Inertial parameter 


and y directions, respectively 

k f 

Absorption coefficient 

x,y 

Cartesian coordinates along the surface and 
normal to it, respectively 
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Greek symbols 

^ Similarity variable 

P Density of the fluid 

l 1 Viscosity of the fluid 

a e Electrical conductivity 

K Thermal conductivity 

15 Kinematic viscosity 

Dimensionless temperature 
® r Temperature ratio parameter 

Superscript 

Derivative with respect to 

Subscripts 

w Properties at the plate 

® Free stream condition 

INTRODUCTION 

The study of boundary layer flow behavior and heat transfer 
characteristics of a Newtonian fluid past a plate embedded 
in a fluid saturated porous medium is important because the 
analysis of such flows finds extensive applications in engi¬ 
neering processes, especially in the enhanced recovery of 
petroleum resources and packed bed reactors (Pal and Shi- 
vakumara [1]). Abetter understanding of convection through 
porous medium can benefit several areas like geophysical 
flow problems, grain storage, catalytic reactor, insulation de¬ 
sign, material handling conveyers, geothermal system, filter¬ 
ing devices etc. In recent years, MHD flow problem have 
become more important industrially. Indeed, MHD bound¬ 
ary layer flow behavior and heat transfer characteristics of 
a Newtonian fluid past a vertical plate embedded in a fluid 
saturated porous medium is significant type of flow having 
considerable practical application in engineering processes. 
The MHD boundary layer flow of an electrically conduct¬ 
ing fluid is encountered in many problems in geophysics, 
astrophysics, engineering applications and other industrial 
areas. Hydromagnetic free convection flows have great sig¬ 
nificance for the applications in the fields of stellar and plan¬ 
etary magnetospheres, aeronautics. Engineers apply MHD 
principle in the design of heat exchangers, pumps, in space 
vehicle propulsion, thermal protection, control and re-entry 
and in creating novel power generating systems. However, 
hydromagnetic flow and heat transfer problems have become 
more important industrially. In many metallurgical processes 
involve the cooling of many continuous strips or filaments 
by drawing them through an electrically conducting fluid 
subject to a magnetic field, the rate of cooling can be con¬ 


trolled. 

Convective heat transfers with thermal radiation are very 
important in the process involving high temperature such as 
gas turbines, nuclear power plant and thermal energy stor¬ 
age etc. in light of these various applications, Hossain and 
Takhar [2] studied the effect of thermal radiation using Ros- 
seland diffusion approximation on mixed convection along 
a vertical plate with uniform free stream velocity and sur¬ 
face temperature. Furthermore Hossain et al. [3-4] studied 
the thermal radiation of a grey fluid which is emitting and 
absorbing radiation in a non-scattering medium. Free con¬ 
vective flows in a saturated porous medium were made by 
Cheng and Minkowycz [5] and Cheng [6]. Wilks [7] dis¬ 
cussed the combined forced and free convection flow along 
a semi-infinite plate extending vertically upwards. Lai and 
Kulacki [8] investigated the coupled heat and mass transfer 
by mixed convection from a vertical plate in a saturated po¬ 
rous medium. Furthermore, some important contribution in 
Darcian and non-Darcian mixed convection about a vertical 
plate were made by Hsu and Cheng [9], Vafai and Tien [10]. 
Soundalgekaret et al. [11] discussed the combined free and 
forced convection flow past a semi-infinite plate with varia¬ 
ble surface temperature. The problem of Darcy-Forchheimer 
mixed convection heat and mass transfer in fluid-saturated 
porous media was studied by Rami et al. [12]. Goren [13] 
was one of the first to study the role of thermophoresis in 
the laminar flow of a viscous and incompressible fluid. Most 
previous studies of the same problem neglected viscous dis¬ 
sipation and thermophoresis. But Gebhart [14] has shown 
that the viscous dissipation effect plays an important role 
in natural convection in various devices that are subjected 
to large variations of gravitational force or that operate at 
high rotational speeds. Analytically study of the non-Darcian 
effects on a vertical plate natural convection in porous me¬ 
dia was made by Hong et al. [15]. Kaviany [16] studied the 
Darcy-Brinkman model to study the effects of boundary and 
inertia forces on forced convection over a fixed impermeable 
heated plate embedded in a porous medium. Chen and Ho 
[17] discussed the effects of flow inertia on vertical, natu¬ 
ral convection in saturated porous media. Kumari et al. [18] 
investigated the non-Darcian effects on forced convection 
heat transfer over a flat plate in a highly porous medium. 
Anghel et al. [19] studied the Dufour and Soret effects on 
free convection boundary layer over a vertical surface em¬ 
bedded in a porous medium. The effects of magnetic field 
and thermal radiation on forced convection flow were made 
by Damesh et al. [20]. Samad and Rahman [21] investigated 
the thermal radiation interaction with unsteady MHD bound¬ 
ary layer flow past a continuous moving vertical porous plate 
immerse in a porous medium with time dependent suction 
and temperature, in presence of magnetic field with radia¬ 
tion. Postelnicu [22] analyzed the influence of magnetic field 
on heat and mass transfer from vertical surfaces in porous 
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media considering Soret and Dufour effects. Srinivasacha- 
rya and Upendar [23] examined the free convection in MHD 
micropolar fluid under the influence of Dufour and Soret 
effects. Duwairi and Damesh [24] studied the natural con¬ 
vection heat and mass transfer by steady laminar boundary 
layer flow over an isothermal vertical flat plate embedded in 
a porous medium. Mukhopadhyay et al. [25] investigated the 
forced convection flow and heat transfer in a porous medium 
using the Darcy-Forchheimer model. 

A new dimension is added to the abovementioned study by 
considering the effects of Darcy-Forchheimer porous media. 
It is well known that Darcy’s law is an empirical formula re¬ 
lating the pressure gradient, the bulk viscous fluid resistance 
and the gravitational force for a forced convective flow in a 
porous medium. Deviations from Darcy’s law occur when 
the Reynolds number based on the pore diameter is within 
the range of 1 to 10 (Ishak et al. [26]). For flow through 
porous medium with high permeability, Brinkman [27] as 
well as Chen et al. [28] argue that the momentum equation 
must reduce to the viscous flow limit and advocate that clas¬ 
sical frictional terms be added in Darcy’s law. Mishra and 
Jena [29] find numerical Solution of MHD flow with vis¬ 
cous dissipation. Vidyasagar et al. [30] discussed heat and 
mass Transfer effects on MHD boundary layer flow over 
a moving vertical porous plate. Mukhopadhyay et al. [31] 
investigated forced convective flow and heat transfer over 
a porous plate in a Darcy-Forchheimer porous medium in 
presence of radiation. Kankanala and Bandari [32] studied 
mixed convection flow of a casson fluid over an exponential¬ 
ly stretching surface with the effects of soret, dufour, ther¬ 
mal radiation and chemical reaction in presence of magnetic 
field. Jena [33] discussed numerical solution of boundary 
layer flow with viscous dissipation in presence of magnetic 
field. Venkateswarlu et al [34] studied effects of chemical 
reaction and heat Generation on boundary layer flow of a 
moving vertical plate with suction and dissipation in pres¬ 
ence of magnetic field. In the present paper, motivated by 
works mentioned above and practical applications, an analy¬ 
sis is made to study the problem of two dimensional MHD 
boundary layer forced convective flow and heat transfer of 
an incompressible fluid past a plate embedded in a Darcy- 
Forchheimer porous. Thermal radiation term is considered 
in the energy equation. Using similarity transformation, the 
governing partial differential equations are transformed into 
a set of self-similar ordinary differential equations, which are 
then solved using Runge-Kutta fourth order method along 
with shooting technique. The numerical results are plotted in 
some figures and the variations in physical characteristics of 
the flow dynamics and heat transfer for several parameters 
involved in the equations are discussed. 


nar boundary-layer flow of a viscous, incompressible, radi¬ 
ating, electrically conducting fluid over a flat plate of very 
small thickness and much larger breadth, embedded in a po¬ 
rous medium. The x-coordinate is measured along the plate 
from its leading edge and the y-coordinate is normal to it. A 
magnetic field B 0 of uniform strength is applied transversely 
to the direction of the flow. In the analysis of flow in porous 
media, the differential equation governing the fluid motion is 
based on Darcy-Forchheimer model, which accounts for the 
drag (represented by the Darcy term) exerted by the porous 
media as well as the inertia effect (represented by the non¬ 
linear Forchheimer term). The magnetic Reynolds number 
is assumed to be small so that the induced magnetic field is 
neglected. The Hall effects, the viscous dissipation and the 
joule heating terms are also neglected. Taking into account 
the thermal radiation term and under the usual boundary lay¬ 
er approximations, the governing equations that are based 
on balance laws of mass, linear momentum and energy (Pai 
[35], Schlichting [36], Bansal [37]) under the influence of 
externally imposed transverse magnetic field (Jeffery [38], 
Bansal [39]) for this investigation can be written as: 


-+- = 0 
Sat dy 


(i) 


Bu B^u v r k f 2 

u^+v-=v — 


dy 


By 


( 2 ) 


u dr . dr _ _1 dq T 

u dx 1 v 8y pCp 8y z pC p 8y 

Accompanied by the boundary conditions: 

y = 0: u = O, v = 0,T = T w 


( 3 ) 


y-t!o:u = u E ,T = T ff (4) 

The governing partial differential equations (l)-(3) can be 
reduced to ordinary differential equations by introducing the 
following transformation: 

n = y (^) ’ 'f = (u„vx>f(n> ( 5 ) 

Where n is the similarity variable and ¥ is the stream func- 
tion defined in the usual way as u = —, v = — — which iden¬ 
tically satisfies (1). 

The transformed non-linear ordinary differential equations 
are: 


f + iff" - kjCf' - 1) - k 2 (f ' “ -1) + Mf' = 0 (6) 


FORMULATION OFTHE PROBLEM ^f 1 + ^ t 1 + ( 0 r -1) 9 ) 3 ] + — (i + ( 9 r - iM 2 (e r - i)e’ + > = o ( 7 ) 

Consider a forced convective, two-dimensional steady lami- 
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Where, 

Kl bb ^ rb x ~ (Porous medium parameter), k 2 = -= (In¬ 
ertial parameter), 

M = (Magnetic parameter), Pr = —^ (Prandtl number), 
N = ———7- (Radiation parameter), 

4pT tE 

T 

0 r = (Temperature ratio of the medium). 

The corresponding boundary conditions are reduced to: 

f(0) = S,f(0) = 0,0(0) = 1 
f (qo) -> 1 , 0(33) -> 0 . ( 8 ) 

DISCUSSIONS 

The system of governing equations (6)-(7) together with the 
boundary condition (8) is non-linear ordinary differential 
equations depending on the various values of the physical 
parameters viz. Porous medium parameter k l9 the Inertial 
parameter k 2 , the Magnetic parameter M, the Prandtl number 
Pr, and the Radiation parameter N. The system of equations 
(6)-(7) is solved by Runge-Kutta fourth order scheme with a 
systematic guessing of f (0) and 0 (0) by the shooting tech¬ 
nique until the boundary conditions at infinity are satisfied. 
The step size Aq = 0.01 is used while obtaining the numerical 
solution and accuracy upto the seventh decimal place i.e. 1 x 
10 4 , which is very sufficient for convergence. The computa¬ 
tions were done by a program which uses a symbolic and 
computer language Matlab. 

RESULTS 

Figure 1 shows the effect of porosity parameter (k x ) on the 
velocity profile. From this plot it is observed that the effect 
of increasing values of porosity parameter is to increases 
the velocity distribution in the flow region. Accordingly, the 
thickness of the velocity boundary layer decreases. In this 
case, horizontal velocity is found to increase with the in¬ 
creasing values of the permeability of the medium. With a 
rise in permeability of the medium, the Darcian body force 
decreases in magnitude (as it is inversely proportional to the 
permeability Darcian resistance acts to decelerate the fluid 
particles in continua. This resistance diminishes as perme¬ 
ability of the medium increases. So progressively less drag 
is experienced by the flow and flow retardation is thereby 
decreased. Hence the velocity of the fluid increases as poros¬ 
ity parameter increases. 

Figure 2 shows the effect of inertial parameter (k 2 ) on veloc¬ 
ity components of fluid velocity. From this plot it is observed 
that the effect of increasing values of inertial parameter is to 
increases the velocity distribution i.e., the thickness of the 


velocity boundary layer decreases. Thus, the non-Dercian 
term (i.e. Forchheimer term) has a very significant effect on 
the velocity profile. 

The impact of the magnetic parameter (M) on the fluid flow 
is very significant in practical point of view. In Figure 3, the 
variation in velocity distribution for several values of mag¬ 
netic parameter is presented. The dimensionless velocity 
f (q) increases with increasing values of magnetic param¬ 
eter. Accordingly, the thickness of the velocity boundary lay¬ 
er decreases. This happens due to the Lorentz force arising 
from the interaction of magnetic and electric fields during 
the motion of the electrically conducting fluid. To reduce ve¬ 
locity boundary layer thickness the generated Lorentz force 
enhances the fluid motion in the boundary layer region. 

In Figures 4 and 5, the effect of suction parameter (S) are 
presented for non-porous (k 1 =0=k 2 ) and porous media (k 1 
=0.1=k 2 ) respectively. We infer from these figures that the 
horizontal velocity increases with increase in suction param¬ 
eter i.e. suction causes to increase the velocity of the fluid. 
Since the effect of suction is to suck away the fluid near the 
wall, the velocity boundary layer is reduced due to suction 
S (S>0). Consequently the velocity increases. Hence the ve¬ 
locity gradient increases with increasing values of suction 
parameter S (S>0). 

Figures 6 and 7 are plotted for the temperature profiles for 
different values of permeability parameter (kj) and inertial 
parameter (k z ) respectively for the given values of N, 0 r , M, 
Pr and S. We observed from both the figures that temperature 
decreases with increase of permeability parameter and iner¬ 
tial parameter within the boundary layer. 

Figure 8, which illustrate the effect of magnetic parameter 
(M) on the temperature profile. We infer from this figure that 
the temperature decreases with an increase in magnetic pa¬ 
rameter. In this case temperature asymptotically approaches 
to zero in the free stream region. 

Figure 9 (in non-porous media) and Figure 10 (in porous me¬ 
dia), which are the graphical representation of the tempera¬ 
ture profiles for different values of the radiation parameter N 
for the given values of the parameters k 1? k 2? 0 r? M, Pr and 
S. We infer from these figures that temperature of the fluid 
decreases with increase in radiation parameter. The increase 
in radiation parameter means the release of heat energy from 
the flow region and so the fluid temperature decreases as the 
thermal boundary layer thickness becomes thinner. 

Figure 11 and 12 depict the effect of suction parameter S 
when k 1 = 0 =k 2 and k 1 =0.1=k z respectively for the given 
values of N, M, 0 r and Pr. We infer from these figures that 
temperature of the fluid decreases with increase in suction 
parameter in both the cases. The thermal boundary layer 
thickness decreases with the suction parameter which causes 
an increase in the rate of heat transfer. The reason of this 
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behavior is that the fluid is brought closer to the surface and 
reduces the thermal boundary layer thickness. 

The effect of temperature ratio parameter 0^ on the tempera¬ 
ture profiles are presented in Figures 13 and 14, in the case 
of non-porous media and porous media respectively. We ob¬ 
served from both the figures that the fluid temperature in¬ 
creases with increase of B r . 9 r >l when T W >T M and this is the 
usual case that is considered in this problem, 9 r <1 repre¬ 
sents the case when T w < T K and this is possible if the porous 
plate is kept on ice pad. Figure 15, which is a representation 
of the local dimensionless coefficient of heat transfer — 0 (0) 

, knows as the Nusselt number for different values of Prandtl 
number (Pr) versus radiation parameter (N) for the given val¬ 
ues of 6 r , M, S. We observed from this figure that the rate of 
heat transfer increases with the value of Pr. 


CONCLUSION 

A mathematical model has been presented for the MHD 
boundary layer forced convective flow and heat transfer of 
an incompressible fluid past a plate embedded in a Darcy- 
Forchheimer porous medium. The governing partial dif¬ 
ferential equations are converted into ordinary differential 
equations by using similarity transformations. The effect of 
several parameters controlling the velocity and temperature 
profiles are shown graphically and discussed briefly. The 
influence of the physical parameters viz. Porous medium pa¬ 
rameter k 1 , the Inertial parameter k 2 , the Magnetic parame¬ 
ter M, the Prandtl number Pr, and the Radiation parameter N 
on dimensionless velocity and temperature profiles were ex¬ 
amined. From the study, following conclusion can be drawn: 

1. The effect of permeability of the medium on a viscous 
incompressible fluid is to increase the fluid velocity by 
reducing the drag on the flow which in turn causes a 
decrease in the temperature field. 

2. It is found that the inertial parameter has a great in¬ 
fluence on decreasing the flow field, whereas its in¬ 
fluence is reversed on the rate of heat transfer. The 
non-Darcian term (i.e. Forchheimer term) has a very 
significant effect on the velocity distribution. 

3. As the magnetic parameter M increases, we can find 
the velocity profile increases in the flow region and 
to decrease the temperature profile. This happens due 
to the Lorentz force arising from the interaction of 
magnetic and electric fields during the motion of the 
electrically conducting fluid. Thus we conclude that 
we can control the velocity field and temperature by 
introducing magnetic field. 

4. The velocity gradient increases with increasing values 
of suction parameter S. Since the effect of suction is to 
suck away the fluid near the wall, the velocity bound¬ 
ary layer is reduced due to suction. Consequently the 
velocity increases. Hence the velocity gradient in¬ 


creases with increasing suction. 

5. The temperature in the boundary layer decreases due 
to suction. Due to thermal radiation, temperature is 
found to decrease. This effect is more pronounced in 
presence of porous media. The combined effects of 
suction and thermal radiation can be used as means of 
cooling. 

6. The fluid temperature increases with increase of tem¬ 
perature ratio parameter 0 p . 

7. The rate of heat transfer increases with the increasing 
values of Prandtl number. This heat transfer is very 
important in production engineering to improve the 
quality of the final product. 
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Figure 2: Velocity profile for various values of k 2 when M=.01, 
N=1, Pr=1, S=.5, 0 r =1.1, ^=0.1. 
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Figure 3: Velocity profile for various values of M when N=1, 
Pr=1, S=.5, Q r =1.1, ^ =0.1, k 2 =0.1 



Figure 4: Velocity profile for various values of S when N=1, 
Pr=1, M=.01, d r =1.1, =0, & 2 =0. 



Figure 5: Velocity profile for various values of S when N=1, 
Pr=1, M=.01, 0 r =1.1, k 2 =0A, k 2 —0.1. 



Figure 7: Temperature profile for various values of & 2 when 
M=.01, N=1, Pr=1, S=.5, k x = 1 . 1 , k x =o.1. 



Figure 8: Temperature profile for various values of M when 
N=1, Pr=1, S=.5, @ r =1.1, ^2=0.1, k 2 - 0.1. 



Figure 9: Temperature profile for various values of N when 
M=01, Pr=1, S=.5, 0 r =1.1, £j = 0, k 2 = o. 



Figure 6: Temperature profile for various values of k x when Figure 10: Temperature profile for various values of N when 
M=.01, N=1, Pr=1, S=.5,^.=1.1, k 2 =QA. M=.01, Pr=1, S=.5, 0 r =1.1, K =0A,k 2 =0.1. 
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Figure 11: Temperature profile for various values of S when 
M=01, Pr=1, N=1,=1.1, £ 1 = 0, ^2 = 0 


Figure 14: Temperature profile for various values of Q when 
M=01, Pr=1, N=1 ,S=.5, ^ = 0.1, & 2 = 0.1 



Figure 12: Temperature profile for various values of S when 
M=.01, Pr=1, N=1, K =1.1, ^i = 0.1, = 0.1 




Figure 15: Variation of — O r'O 'i with the Radiation parameter 
N for various values of Pr when M=.01, S=.5, ^ =1.1, = 0 , 

^2 = 0 . 


Figure 13: Temperature profile for various values of k x when 
M=.01, Pr=1, N=1 ,S=.5, k x = 0, 6 r = 0. 


1 


Int J Cur Res Rev | Vol 9 • Issue 11 • June 2017 
































